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ABSTRACT 

We consider the luminosity and environmental dependence of structural parameters of 
lenticular galaxies in the near-infrared K band. Using a two-dimensional galaxy image de- 
composition technique, we extract bulge and disk structural parameters for a sample of 36 
lenticular galaxies observed by us in the K band. By combining data from the literature for 
field and cluster lenticulars with our data, we study correlations between parameters that char- 
acterise the bulge and the disk as a function of luminosity and environment. We find that 
scaling relations such as the Kormendy relation, photometric plane and other correlations in- 
volving bulge and disk parameters show a luminosity dependence. This dependence can be ex- 
plained in terms of galaxy formation models in which faint lenticulars (Mt > —24.5) formed 
via secular formation processes that likely formed the pseudobulges of late-type disk galax- 
ies, while brighter lenticulars (Mr < —24.5) formed through a different formation mecha- 
nism most likely involving major mergers. On probing variations in lenticular properties as 
a function of environment, we find that faint cluster lenticulars show systematic differences 
with respect to faint field lenticulars. These differences support the idea that the bulge and 
disk components fade after the galaxy falls into a cluster, while simultaneously undergoing a 
transformation from spiral to lenticular morphologies. 
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1 INTRODUCTION 

Lenticular (SO) galaxies were originally conceived as a morpho- 
logical transition class between ellipticals and early-type spirals by 
Hubble (1936). These galaxies have disks spanning a wide range 
in luminosity, contributing between five to fifty percent of the to- 
tal galaxy light. They are distinguished from spiral galaxies by the 
lack of conspicuous spiral arms. For these reasons, lenticulars can 
be thought of as a population which is intermediate between ellip- 
ticals and early type spirals. Indeed, the placement of lenticulars by 
Hubble (1936) on the tuning fork diagram clearly implies such a 
relationship. In several observable properties such as bulge-to-disk 
luminosity ratio, star formation rate and colour, lenticulars are, on 
average, intermediate between ellipticals and early type spirals. 

Detailed study of individual lenticular galaxies indicates that 
the situation is more complex in reality. It has been suggested (van 
den Bergh 1994) that there are different, but overlapping, sub- 
populations amongst the lenticulars. Exploration of formation sce- 
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narios for lenticulars using theory and numerical simulations also 
suggest that lenticulars may have formed in different ways. They 
could be of primordial origin forming rather rapidly at early epochs, 
or could have been formed by the slow stripping of gas from spirals, 
which changes the morphology (Abadi, Moore & Bower 1999), or 
through the mergers of unequal-mass spirals (Bekki 1998). The two 
main components of lenticulars - the bulge and the disk - may have 
their own different and possibly independent formation history. Ob- 
servationally, the disks seem to be younger than the bulges in both 
spiral and lenticular galaxies (Peletier & Balcells 1996), with the 
age difference between the two components larger in lenticulars 
(Bothun & Gregg 1990). Bars in lenticular galaxies are also of in- 
terest in many studies as they provide a clue about the evolutionary 
history of these galaxies. Recent studies reveal that bars in lenticu- 
lar galaxies are shorter, less massive and have smaller bar torques 
(Buta et al. 2006; Laurikainen et al. 2005; Laurikainen et al. 2006; 
Gadotti, D. A. et al. 2007). 

A detailed multiband study of the morphology of representa- 
tive samples of lenticulars in different environments, and compar- 
ison of their properties with those of ellipticals, and with bulges 
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Table 1. Best fit Bulge and Disk parameters for our sample. 
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UGC 00080 


0.01019 


-3 


11.54 


0.45 


1.63 


0.001 


15.69 


2.33 


0.505 


0.19 


-24.99 


UGC 00491 


0.01664 


-1 


10.38 


4.95 


3.82 


0.257 


17.65 


5.08 


0.231 


0.65 


-25.53 


UGC 00859 


0.00711 





11.34 


1.39 


3.36 


0.141 


16.40 


1.22 


0.403 


0.43 


-23.25 


UGC 00926 


0.01467 


-2.5 


09.74 


9.81 


4.39 


0.348 


14.93 


0.80 


0.353 


0.93 


-26.17 


UGC 01250 


0.01235 


-2 


09.51 


1.92 


3.79 


0.397 


15.92 


2.66 


0.696 


0.32 


-25.24 


UGC 01823 


0.01332 


-3 


11.41 


1.96 


2.28 


0.416 


16.81 


6.56 


0.398 


0.34 


-26.34 


UGC 01964 


0.01732 


-2 


10.06 


4.03 


3.96 


0.103 


15.60 


1.86 


0.585 


0.58 


-25.33 


UGC 02039 


0.01505 


-2 


10.18 


1.92 


3.37 


0.007 


17.19 


5.52 


0.644 


0.29 


-25.52 


UGC 02187 


0.01605 


-2 


13.11 


0.45 


0.76 


0.276 


15.13 


1.99 


0.554 


0.14 


-25.17 


UGC 02322 


0.01447 


-1 


10.51 


4.18 


3.86 


0.008 


16.84 


2.08 


0.015 


0.75 


-24.89 


UGC 03 178 


0.01578 


-2 


09.84 


2.55 


3.89 


0.014 


14.76 


1.17 


0.479 


0.47 


-24.90 


UGC 03452 


0.01876 


-2 


10.53 


3.31 


3.47 


0.242 


15.91 


2.89 


0.516 


0.46 


-25.70 


UGC 03536 


0.01564 


-2 


13.18 


0.55 


0.99 


0.012 


14.80 


2.03 


0.576 


0.10 


-25.48 


UGC 03567 


0.02016 


-2 


09.96 


4.87 


4.15 


0.062 


15.99 


1.06 


0.503 


0.87 


-25.04 


UGC 03642 


0.01500 


-2 


11.94 


1.03 


2.02 


0.010 


17.74 


8.45 


0.206 


0.17 


-25.72 


UGC 03683 


0.01908 


-2 


12.13 


1.56 


1.83 


0.283 


17.01 


5.69 


0.171 


0.37 


-25.91 


UGC 03699 


0.01947 


-2 


10.35 


3.13 


3.32 


0.336 


16.54 


1.65 


0.489 


0.87 


-25.36 


UGC 03792 


0.01908 





10.17 


4.41 


3.74 


0.224 


18.15 


5.71 


0.225 


0.72 


-25.66 


UGC 03824 


0.01787 


-2 


11.40 


1.52 


2.71 


0.003 


17.43 


3.18 


0.002 


0.52 


-24.52 


UGC 04347 


0.01494 


-2 


10.68 


4.33 


3.34 


0.243 


16.34 


1.01 


0.528 


0.95 


-25.57 


UGC 04767 


0.02413 


-2 


09.05 


7.65 


4.74 


0.105 


14.98 


0.27 


0.001 


0.99 


-25.55 


UGC 04901 


0.02810 


-2 


13.30 


3.06 


1.78 


0.305 


18.08 


14.28 


0.043 


0.27 


-26.70 


UGC 05292 


0.00511 


-2 


09.22 


0.26 


3.08 


0.008 


17.28 


2.01 


0.282 


0.21 


-23.07 


UGC 06013 


0.02190 


-2 


11.29 


5.10 


3.68 


0.325 


15.41 


1.00 


0.001 


0.78 


-24.89 


UGC 06389 


0.00663 


-2 


11.71 


3.00 


3.36 


0.378 


13.88 


0.21 


0.334 


0.89 


-23.75 


UGC 06899 


0.02249 


-2 


09.79 


5.11 


4.06 


0.060 


16.71 


0.53 


0.001 


0.99 


-25.35 


UGC 07142 


0.00328 


-2 


10.48 


0.35 


2.22 


0.219 


15.96 


1.32 


0.503 


0.33 


-23.66 


UGC 07473 


0.00414 


-2 


11.76 


0.24 


1.20 


0.239 


14.55 


1.20 


0.725 


0.12 


-24.57 


UGC 07880 


0.00388 


-3 


11.50 


0.68 


2.39 


0.330 


14.29 


0.63 


0.747 


0.34 


-23.74 


UGC 08675 


0.00361 


-2 


09.18 


0.25 


3.35 


0.020 


17.14 


1.48 


0.011 


0.19 


-22.53 


UGC 09200 


0.01081 


-2 


12.33 


1.35 


2.23 


0.117 


17.83 


3.10 


0.001 


0.57 


-24.15 


UGC 09592 


0.01791 


-2 


08.69 


1.62 


3.62 


0.315 


16.30 


2.33 


0.123 


0.65 


-25.30 


UGC 11356 


0.00820 


-2.5 


11.64 


0.47 


1.64 


0.080 


15.16 


1.46 


0.020 


0.24 


-24.57 


UGC 11972 


0.01463 


-2.5 


10.29 


8.62 


4.08 


0.441 


14.24 


0.45 


0.002 


0.95 


-26.00 


UGC 12443 


0.01463 


-2 


10.83 


3.34 


3.44 


0.097 


19.36 


0.04 


0.002 


1.00 


-24.61 


UGC 12655 


0.01727 


-2 


10.34 


1.14 


2.87 


0.191 


16.86 


4.77 


0.400 


0.24 


-25.47 



Notes. Column (1) gives the UGC catalogue number, column (2) and column (3) give redshift and morphological type respectively from NED, 
columns (4), (5), (6) and (7) give unconvolved bulge central surface brightness, bulge effective radius, Sersic index and bulge ellipticity respectively, 
columns (8) (9) and (10) give unconvolved disk central surface brightness, disk scale length and disk ellipticity respectively, column (11) gives the 
bulge-to-total luminosity ratio and column (12) gives the absolute magnitude in the K band. 



and disks of spirals will be important in addressing these issues 
observationally. However, comparing the predictions of models to 
observations is complicated for two main reasons: (1) The mod- 
els often use simplifying assumptions that may not hold for real 
galaxies and (2) Models often do not make firm predictions about 
directly observable quantities. Nevertheless, the statistical proper- 
ties of galaxy ensembles can be compared to model predictions. 
Such a quantitative comparison between bulge and disk properties 
of galaxies is considerably simplified if one assumes simple ana- 
lytic profiles for the light distribution of the bulge and the disk. In 
practice, the Sersic r 1 /™ law (Sersic 1968) adequately represents 
the bulge surface brightness distribution while an exponential best 
represents the disk surface brightness distribution. The total galaxy 
light is simply the sum of a Sersic bulge and an exponential disk. 

In this work, we obtain and report structural parameters for 
36 lenticular galaxies observed in the near-infrared K band. Us- 



ing results from the bulge-disk decomposition, we examine corre- 
lations among the bulge and disk parameters and discuss implica- 
tions to models of lenticular formation. Our goal is mainly to study 
the constraints placed on bulge formation in lenticular galaxies by 
these parameter correlations. Specifically, we want to investigate, 
in greater detail, the discovery of Barway et al. (2007) who found 
two distinct populations of bulges in lenticulars that were fainter 
and brighter than a threshold luminosity. 

This paper is organised as follows: the sample, details of near- 
infrared observations, data reduction technique and photometric 
calibration are summarised in § 2. The two-dimensional decompo- 
sition analysis is described in § 3. In § 4, we discuss the luminosity 
and environment dependence of lenticular galaxy properties. 

Throughout this paper, we use the standard concordance cos- 
mology with Q.m = 0.3, S7a = 0.7, and hum = 0.7. 
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2 OBSERVATIONS AND DATA REDUCTION 
2.1 Our sample 

Our sample consists of a set of 36 bright field lenticular galax- 
ies from Barway et al. (2005). The original sample contained 40 
lenticulars galaxies, selected from the Uppsala General Catalogue 
(UGC), with apparent blue magnitude brighter than ms = 14, di- 
ameter D25 < 3 arcmin and declination in the range 5 < 8 < 64° . 
The sample, while not complete, is representative of bright field 
lenticulars. 

We obtained images of the sample galaxies in the near in- 
frared K band with the Observatorio Astronomico National 2.1-m 
telescope at San Pedro Martir, Mexico. The CAMILA instrument 
(Cruz-Gonzalez et al. 1994), which hosts a NICMOS 3 detector 
of 256x256 pixel format, was used in the imaging mode. Each K 
band observing sequence consisted of 10 exposures, six on the ob- 
ject and four on the sky. The net exposure times were, typically, 10 
minutes per galaxy. A series of twilight and night-sky images were 
taken for flat-fielding purposes. The K observations were carried 
out in four runs in December 2000, March 2001, October 2001 and 
March 2002. The data reduction procedure for the K images in- 
volved subtraction of the bias and sky frames, division by flat field 
frames, registration of the images to a common co-ordinate system 
and then stacking of all the images of a given galaxy. All image re- 
ductions were carried out using the Image Reduction and Analysis 
Facility (IRAF 1 ) and the Space Telescope Science Data Analysis 
System (STSDAS 2 ). Standard fields were observed in order to en- 
able accurate photometric calibration of our K band observations. 
Reddening corrections due to Galactic extinction and K-correction 
were applied to individual galaxies. Full details on sample selec- 
tion, observation, and data reduction procedures can be found in 
Barway et al. (2005). 



2.2 Comparison samples 

We supplement our sample with data from Bedregal et al. 2006 
(hereafter BAM06) and additional data provided by Bedregal & 
Aragon-Salamanca in electronic form (private communication) 
who used the Two Micron All Sky Survey (2MASS; Jarrett et 
al. 2003) data for a structural analysis of a sample of 49 lenticu- 
lar galaxies. These are relatively faint objects with sufficient rota- 
tional support for the disks. Addition of these data complements 
our sample in two ways: (a) it provides a low luminosity exten- 
sion to the galaxies in our sample and (b) it provides lenticulars in 
different environments i. e. , galaxies from the Coma (14 galaxies), 
Virgo (8 galaxies), and Fornax (6 galaxies) clusters along with 21 
field lenticulars. Whenever a comparison with other morphological 
classes is appropriate, we have used results from analyses of the 
following samples: 

• A sample of 42 elliptical galaxies from the Coma cluster ob- 
served in K band by Mobasher et al. (1999) as analysed by Khos- 
roshahi et al. (2000b). 

• A sample of 26 early type spiral galaxies in the field observed 
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by Peletier & Balcells (1997) in the K band as analysed by Khos- 
roshahi et al. (2000b). 

• A sample of 40 late-type spiral galaxies observed and analysed 
by Mollenhoff & Heidt (2001) in the K band. 

All the above analyses model the galaxy light using the Sersic 
function for the bulge and an exponential function for the disk. 
In all cases, the bulge-disk decomposition is performed using the 
full 2D image of the galaxy. Incidentally, for all samples except 
that of BAM06 and Mollenhoff & Heidt (2001), the decomposition 
has been performed with the same code fitgal (Wadadekar et al. 
1999). Details of the sample selection, observation, data reduction 
and bulge-disk decomposition of these samples may be found in 
the references cited. 



3 ANALYSIS 

3.1 2-d image decomposition 

Extracting the structural parameters of a galaxy requires the sepa- 
ration of the observed light distribution into bulge and disk compo- 
nents. There is considerable variation in the details of the decompo- 
sition techniques proposed by various researchers. In recent years, 
methods that employ two dimensional fits to broad band galaxy 
images have become popular (e. g. Wadadekar, Robbason & Kem- 
bhavi 1999; Peng et al. 2002; Simard et al. 2002; de Souza et 
al. 2004). Most of these decomposition techniques assume specific 
surface brightness distributions like the Sersic law for the bulge and 
an exponential distribution for the disk. 

Our decomposition procedure is a full two dimensional 
method that uses information from all pixels in the image, fitgal es- 
sentially involves a numerical solution to a signal-to-noise ( here- 
after S/N) weighted \ 2 minimisation problem. We achieve this 
minimisation using the Davidon-Fletcher-Powell variable metric 
algorithm included as part of MINUIT - a multidimensional min- 
imisation package from CERN. Our technique involves iteratively 
building two dimensional image models that best fit the observed 
galaxy images, with the quality of the fit quantified by the \ 2 value. 
We compute weights for the x 2 function using the S/N ratio at 
each pixel of the galaxy image. We convolve the model image with 
the measured point spread function from the galaxy frame before 
the \ 2 ' s computed. Details of the accuracy and reliability of the 
decomposition procedure as assessed by simulations, are provided 
in Wadadekar et al. (1999). 

At near-infrared wavelengths, dust related absorption and 
emission from patchy star forming regions are both weak relative 
to optical wavelengths. The smooth, featureless light profiles of 
galaxies in the near-infrared are very convenient for modelling us- 
ing simple analytic functions to represent the galaxy components. 
Use of K band data is especially important for our sample, as it 
includes several dusty galaxies. In our bulge-disk decomposition, 
we have the following quantities as free parameters: (1) lb(0): the 
central bulge intensity, in counts, which can later be converted to 
magarcsec -2 using the photometric calibration (2) r e : the half 
light radius of the bulge in pixels (3) e b . the ellipticity of the bulge 
(4) n: the bulge Sersic index (5) Jd(0): the central intensity of the 
disk in counts (6) r^: the scale length of the disk in pixels (7) e^: 
the ellipticity of the disk 

With these definitions, the Sersic bulge intensity distribution 
can be written as 

hui ge (x,y) = l b (0) e - 2 - 303b ^ r >>^ /r ° )1/n , (1) 
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Figure 1. Distribution of total absolute magnitude (My) in K band for our 
sample (solid line) and for BAM06 lenticulars (dotted line). The vertical 
dashed line corresponds to total absolute magnitude My = -24.5, which we 
use to divide low- and high-luminosity lenticulars. 



rtuige = \Jx 2 + y 2 /(l - e 6 ) 2 , 

where x and y are the distances from the centre of the galaxy along 
the major and minor axis respectively and b n is a function of n and 
the root of an equation involving the incomplete gamma function. 
However, following Khosroshahi et al. (2000b), b n can be approx- 
imated as a linear function of n, accurate to better than one part in 
10 5 , by 

b n = 0.868242 n- 0.142058. 

For n — 4, which corresponds to de Vaucouleurs law, 64 = 3.33. 

The projected disk profile is represented by an exponential dis- 
tribution, 



Idisk(x,y) = ld(0)e 



' r disk/ r d 



(2) 



rdxsk = V x ' 2 + y 2 /(! 



ed 



The ellipticity of the disk in the image is due to projection effects 
alone and is given by 



■ cos(i), 



(3) 



where i is the angle of inclination between the line of sight and the 
normal to the disk plane. 

The bulge-to-disk luminosity ratio is a dimensionless parame- 
ter which is commonly used as a quantitative measure for morpho- 
logical classification of galaxies. For a Sersic bulge and an expo- 
nential disk it is given by 



(Rim nr(2n) fh(0)\ fr e 

{B/V)n - (2.303M 2 " U(0)J U 



(4) 



We use the bulge-to-total luminosity ratio B/T = (B/D)/{1 + 
(B/D)) in this paper, since it spans a restricted range of to 1. 



3.2 Extraction of structural parameters using 2-d image 
decomposition 

Although the K band observations are advantageous for their rela- 
tive lack of absorption related problems, the high sky background 
in K band poses some difficulty in extracting the global photomet- 
ric parameters accurately, as it limits the region where S/N is high 
enough for the pixels to be usable. Even a small error in the sky 
estimation can lead to spurious results because of the relatively low 
S/N in the images. For example, if the sky is under-estimated, say 
by 5%, the residual background can lead to a spurious detection 
of a large circular disk with a non-zero 7d(0). Even when a 'bulge 
only' model is used, a wrong assessment of the background can 
lead to serious mis-estimation of the morphological parameters. 

In our K band images, the sky background is automatically 
subtracted from the many individual sky frames for each galaxy 
during the preprocessing stage (see Barway et al. 2005 for details). 
However, if sky is removed, noise statistics calculations at low S/N 
(where the sky dominates) are grossly inaccurate. To avoid this 
problem, we added back the sky background that had been sub- 
tracted from each galaxy image. 

Our bulge-disk decomposition code fitgal (Wadadekar et al. 
1999) allows the estimation of the background simultaneously with 
the estimation of morphological parameters. However, such a mea- 
surement can lead to unstable values for the parameters, as the 
background counts are much higher than the signal even in the cen- 
tral region of the galaxies. To avoid this instability, we had a run 
for each galaxy in which the aim was to estimate the background 
rather than to obtain the morphological parameters accurately. To 
this end, we first ran the code on a selected region of each galaxy 
image, dominated by the sky background. We compared the value 
of the background returned by the fitting procedure, with that es- 
timated from regions of the image frames free of the galaxy and 
from the observed sky images to ensure that these values are all 
consistent with each other. We then adopted the sky value returned 
by the code as the best representation of the sky background, since 
this uses information from all the relevant pixels of the image, and 
fixed the background parameter to this value during all subsequent 
runs of fitgal. 

We fitted ellipses to the isophotes of each galaxy of our sam- 
ple, and obtained the one dimensional surface brightness profile 
along the major axis, using the STSDAS task ellipse. For this task, 
sky subtracted galaxy frames were used. The fits provided us with 
model independent surface brightness distributions, which could 
be compared to the surface brightness distribution predicted by our 
best fit bulge plus disk models. In order to obtain good initial val- 
ues for our two dimensional fits, we fitted the inner part of one 
dimensional profiles obtained above with de Vaucouleurs law pro- 
files (taking care to avoid the region affected by the PSF) and the 
outer part with the exponential law. This provided us with approx- 
imate values of lb(0) and r c for bulge and /d(0) and r<j for disk, 
which we used as input to obtain a full two dimensional Sersic and 
exponential model fit to each galaxy in the sample. During the fit 
for each galaxy, we fixed the background to the value obtained 
as described above. Of the 40 galaxies in our sample, we could 
not get satisfactory fits for three galaxies, UGC 03087, UGC 1 1 178 
and UGC 11781. UGC 03087 is a strong radio source with an op- 
tical jet visible in the image make the fit unreliable. For galaxies 
UGC 11178 and UGC 11781, our images suffer from poor S/N 
ratio and we are not able to get satisfactory fits. UGC 07933 is clas- 
sified as an elliptical in the RC3 catalogue. We therefore exclude 
these four galaxies from subsequent discussions, leaving us with a 
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sample of 36 galaxies. We have listed in Table 1 the best fit bulge 
and disk parameters for all the 36 lenticular galaxies in the sample. 

BAM06 used images in the K band from the 2MASS sur- 
vey to obtain bulge and disk parameters for the 49 galaxies in 
their sample, using the gim2d decomposition code by Simard et 
al. (2002), assuming Sersic and exponential laws for the bulge and 
disk light distributions, as we have done with our sample. Barway 
et al. (2007) have shown (see their Figure 2a) that there are three 
clear outliers among the cluster lenticulars of the BAM06 sample. 
Inspection of the 2MASS K band images of the corresponding 
galaxies, shows that one of these (ES0 358-G59) has poor S/N 
ratio, while the other two (NGC4638 and NGC4787) are obvi- 
ously disk-dominated systems, which are likely to have disk scale 
lengths larger than those reported by BAM06. Following Barway 
et al. (2007), we omit these three outliers from further discussion, 
while noting that our conclusions are not significantly changed by 
this omission. 

It must be noted that the depth of our K band images is con- 
siderably more than the depth of 2MASS images used by BAM06. 
We have a typical exposure time of 10 minutes per galaxy as op- 
posed to 7.8 sec/frame in 2MASS. This combined with the fact that 
we used a somewhat larger telescope for obtaining our data allows 
us to reach a depth of typically 22 mag arcsec~ 2 , corresponding 
to an error in photometry of 0.1 mag arcsec~ 2 . The 2MASS data 
used by BAM06 reach a typical depth of 20 mag arcsec~ 2 for the 
same error. This implies that our images are about 2 magnitudes 
deeper than the 2MASS data. Given the significant difference in 
typical depth, bulge disk decompositions of the low S/N 2MASS 
data may, in principle, be systematically affected. To probe the ex- 
tent of this effect, we obtained 2MASS images for all the galaxies 
in our sample and performed the decomposition using our tech- 
nique on these data. For 80% of galaxies, we found a reasonable 
match between parameters extracted from our data and the 2MASS 
data. For 20% of galaxies our fitting procedure failed when using 
2MASS data, indicating that the low S/N was affecting our ability 
to extract parameters reliably for these galaxies. However, such an 
inability to obtain reliable fits is unlikely to affect the parameter 
extraction procedure of BAM06, which is based on the Metropolis 
algorithm (Metropolis et al. 1953), as implemented in the gim2d 
software (Simard et al. 2002). gim2d is a robust, albeit slow, bulge- 
disk decomposition software that obtains the global minimum of a 
fit, in almost any situation. It is thus well suited for parameter ex- 
traction for relatively shallow data like that of BAM06, where the 
low S/N is likely to make fitting difficult. 

Another possible concern is whether a linear combination of a 
bulge and a disk is an adequate mathematical formulation to model 
the light distribution of galaxies in our sample e.g. if a significant 
non-axisymmetric bar (or similar structure) were present for the 
galaxies in our sample, that would affect the parameters that we 
extract (see Laurikainen et al. 2005 for a discussion on the fitting of 
non-axisymmetric components). It is fortunately easy to detect the 
presence of non-axisymmetric structure by examining the residuals 
of our fitting, which are computed as the difference between the 
galaxy image and best fit model. Since the analytic functions we 
use to model the bulge and disk are axisymmetric, the bars are not 
modeled in our scheme. They remain in the residual image, and can 
be visually seen if they are bright enough. For the 36 galaxies in our 
sample, the RC3 catalog indicates that 5 may have bars. We find 
from the residuals that only one galaxy (UGC 80) has a discernable 
bar. For the BAM06 galaxies, the RC3 catalog indicates that 14 out 
of 47 galaxies may have bars. We find non-axisymmetric structures 
in 1 1 of these 15 galaxies. For these galaxies, extracted parameters 
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Figure 2. Dependence of the mean surface brightness within bulge effective 
radius ((/ib(< fc))) on bulge effective radius r c (a) for bright (as open 
circles) and faint lenticulars (filled circles). Coma ellipticals (as dots) are 
overplotted for comparison, (b) for bright lenticulars with bars (as crosses) 
and faint lenticulars with bars (as dashed) as classified in RC3 catalogue. 
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Figure 3. An edge-on view of the photometric plane (PP) for lentic- 
ulars. Top: PP for bright lenticulars; log n = 0.251 ±0.045 log r c - 
0.094±0.012^t b (0) + 1.377±0.138. Bottom: PP for faint lenticulars; logn 
= 0.185±0.029 log r c - 0.041±0.004 ,u b (0) + 0.935±0.046. 



are likely to be somewhat inaccurate. For the analysis subsequently 
reported in this paper, such inaccuracy is of concern only if it leads 
to systematic offsets between barred and non-barred galaxies. We 
demonstrate in Section 4.1.1 that such obvious systematic offsets 
do not exist, implying that inaccurate parameter estimation for a 
small fraction of galaxies will not alter the results of this work. 
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Table 2. photometric plane coefficients for bright & faint lenticulars. 



Lenticular Sample 


a 


b 


c 


N 


rmsn 


rms 


Bright 


0.251 ±0.045 


-0.094±0.012 


1.377±0.138 


38 


0.037 


0.036 


Faint 


0.185±0.029 


-0.041 ±0.004 


0.935±0.046 


44 


0.039 


0.038 


Bright with restricted n 


0.150±0.0213 


-0.056±0.007 


1.050±0.067 


29 


0.019 


0.019 



Notes, a and b are the coefficients of logr e and /ib(0) respectively, c is the constant, N is the number of 
galaxies, rms n is the r.m.s. scatter measured along the log n axis and rms is the r.m.s. scatter measured 
in a direction orthogonal to the best-fit photometric plane. 



4 CORRELATIONS AND DISCUSSION 
4.1 Luminosity dependence 

Barway et al. (2007) have presented evidence to support the 
view that the formation history of lenticular galaxies depends upon 
their luminosity. According to this view, low-luminosity lenticular 
galaxies likely formed by the stripping of gas from the disk of late- 
type spiral galaxies, which in turn formed their bulges through sec- 
ular evolution processes. On the other hand, more luminous lenticu- 
lars likely formed at early epochs through a rapid collapse followed 
by rapid star formation. 

As we mentioned earlier, the BAM06 galaxies complement 
our sample in two ways: they extend to fainter luminosities and pro- 
vide lenticulars in different environments. In Figure 1 we show the 
distribution of total absolute magnitude (Mr) in the K band for the 
combined sample, which is seen to span a wide range in luminosity. 
We divide the combined sample into faint and bright groups, using 
Mt = —24.5 as a boundary. The bright sample has 37 lenticulars 
while 46 lenticulars belong to the faint sample according to this 
luminosity division. The boundary at Mt = —24.5 is somewhat 
arbitrary but our results do not critically depend on small shifts in 
the dividing luminosity. For instance, changing this value by half a 
magnitude on either side maintains correlations significant at least 
at the 95% level. 

Using the luminosity division at Mt = —24.5, Barway et al. 
(2007) reported markedly different correlations between bulge ef- 
fective radius (r ) and disk scale length (r<j) for bright and faint 
lenticulars. A positive correlation of bulge and disk sizes is ex- 
pected if the bulge grows over time through secular evolution pro- 
cesses (see review by Kormendy & Kennicutt 2004). No such corre- 
lation is expected if the bulge formed via merger related processes. 
Barway et al. (2007) found this positive correlation between bulge 
and disk sizes of low luminosity lenticulars; such a correlation was 
not observed for bright lenticulars. Moreover, they found that the 
correlation holds for faint galaxies, irrespective of whether they are 
situated in a field or cluster environment. 

If the differences between low and high luminosity lenticu- 
lars are indeed fundamental, then there should be systematic dif- 
ferences seen between the two populations in other correlations 
among bulge and disk parameters. We study such correlations in 
the following sections. 



4.1.1 Kormendy relation and the effect of bars 

For samples of large elliptical galaxies, the central surface 
brightness pib(0) (or, equivalently, the mean surface brightness 
((Mb(< r c )}) within r c ) is correlated with log r c . This was first 
reported by Kormendy (1977) and is known as the Kormendy re- 
lation. In Figure 2(a) we show the Kormendy relation for bright 
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Figure 4. Photometric Plane for bright lenticulars (open circles) with re- 
stricted n value. Faint lenticulars (filled circles) are over plotted on the same 
photometric plane. These show a higher scatter than the bright ones. 

(as open circles) and faint lenticulars (as filled circles). The bright 
lenticulars show tight correlation with Pearson correlation coeffi- 
cient r — 0.95 with significance greater than 99.99 %. The Ko- 
rmendy relation for faint lenticulars shows considerably greater 
scatter although the Pearson correlation coefficient is 0.75 with 
significance greater than 99.99 %. Figure 2(b) show the Kormendy 
relation for bright (as crosses) and faint (as dashes) lenticulars clas- 
sified as barred lenticulars in RC3 catalogue to probe for any sys- 
tematic effect caused due to presence of bars. It is clear from the 
figure that barred and non-barred show no systematic offset. Only 
three galaxies are classified as barred in our bright galaxy sample 
(indicated by crosses). Even after removing these three galaxies, 
the correlation is nearly unchanged with Pearson correlation coef- 
ficient r = 0.96 with significance greater than 99.99 %. In the faint 
lenticular sample, there are as many as 16 galaxies classified as 
barred. Nevertheless, neglecting these 16 galaxies does not affect 
the correlation which has r — 0.74 with significance greater than 
99.99 %. For all other correlations subsequently reported in this 
paper, we find that excluding barred galaxies does not significantly 
alter the correlation coefficient or its significance. Hence, in further 
discussion, we do not differentiate between barred and non-barred 
lenticulars, and include both types in the analysis. 

The Kormendy relation is known to be a projection of the fun- 
damental plane (Djorgovski & Davis 1987; Dressier et al. 1987) 
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of galaxies. A tight Kormendy relation indicates (near) virialized 
bulges like those found in elliptical galaxies. In Figure 2a, ellipti- 
cals from the Coma cluster are overplotted as dots. As expected, 
these ellipticals show a tight correlation with Pearson correlation 
coefficient r = 0.81 with significance greater than 99.99 %. The 
slopes of the best-fit Kormendy relation for bright lenticulars and 
Coma ellipticals are similar. This similarity coupled with the sys- 
tematic difference in scatter in the Kormendy relation between 
bright and faint lenticulars is an indicator of the more virialized 
state of the bright lenticular bulges and their close relation to ellip- 
ticals. 



4.1.2 Photometric Plane 

Analogous to the fundamental plane (FP), ellipticals and the bulges 
of early type spiral galaxies obey a single planar relation of the form 
log n — a log r c + 6 /i b (0)+c which Khosroshahi etal. (2000a,b) 
called the photometric plane (PP). Assuming that galaxy bulges 
behave as spherical, isotropic, one-component systems, there is a 
unique specific entropy that can be associated with every galaxy 
(Lima-Neto et al. 1999). This in turn, implies a relation between 
bulge parameters that is equivalent to the observed photometric 
plane. An edge-on view of the best-fit PP is shown in the top and 
bottom panels of Figure 3 for bright and faint lenticulars respec- 
tively. Table 2 lists the coefficients for best-fit PP relation. It is im- 
mediately evident from Figure 3 that the PP for faint lenticulars 
spans a limited range of n. We restrict our bright lenticular sample 
to the range spanned by the faint lenticulars and then obtain the best 
fit PP for this restricted sample. The rms scatter for the (restricted) 
bright galaxies is considerably smaller than that for the faint galax- 
ies (see Table 2). In Figure 4 we show the PP for these two sets of 
galaxies. Visually, the higher scatter of the faint galaxies about the 
PP is obvious. The PP is a natural state for galaxies that formed like 
ellipticals (at early epochs in a burst of rapid star formation, with 
stellar orbits that are well relaxed). The fact that bulges of bright 
lenticulars have a tight PP is an indicator of homology with ellipti- 
cals. Faint lenticulars seem to belong to a separate class. 



4.1.3 Correlations involving the Sersic index (n) 

The restricted range spanned by the Sersic index in faint lenticular 
galaxies also affects other correlations involving this index. Fig- 
ure 5 is a plot of the Sersic index (n) as a function of effective radius 
(r e ) which again shows systematic differences between bright and 
faint lenticulars. Bright lenticulars are well correlated having Pear- 
son correlation coefficient r — 0.79 with significance greater than 
99.99 % but faint lenticulars do not show a significant correlation. 
A strong correlation exists between n and the bulge central surface 
brightness (/ib(0)), as shown in Figure 6 for the brighter lentic- 
ulars. The Pearson correlation coefficient in this case is 0.84 at a 
significance level better than 99.99%. Faint lenticulars also exhibit 
good correlation having Pearson correlation coefficient r = 0.58 at 
significance 99.96 %, but have a markedly different slope. Compar- 
ing Figure 5 and Figure 6 with Figure 3 and Figure 4 respectively 
of Khosroshahi et al. (2000b), we observe a close correspondence 
in the distribution of the bulges of bright lenticulars and early type 
spirals. Bulges of early type spirals, in turn, have a close corre- 
spondence with ellipticals (Khosroshahi et al. 2000b). On the other 
hand, faint lenticulars seem to be a different population. 

The bulge-to-total luminosity ratio (B/T) increases with 
Sersic index n as shown in Figure 7 for brighter lenticulars (r = 
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Figure 5. Sersic index n plotted against bulge effective radius r e in kpc. 
Symbols are as in Figure 2. 

0.81 with significance greater than 99.99 %). Faint lenticulars show 
a weak correlation between log (B/T) and log n. Again, the 
trends for bright lenticulars are consistent with those found by 
Khosroshahi et al. (2000b) in early-type spirals. Bright lenticu- 
lars plotted in Figure 3 of Barway et al. (2007) include 5 galaxies 
that are obvious outliers in the anti-correlation of bulge and disk 
sizes. These bright lenticulars viz. UGC 80, UGC 2187, UGC 3536, 
UGC7473 and UGC 11356 do follow the correlations involving 
the Sersic index. However, they are all characterised by a low 
Sersic index, low bulge central surface brightness, low bulge ef- 
fective radius and low bulge-to-total luminosity ratio. This places 
these lenticulars in the bottom left quadrant of Figures 5, 6 and 
7. Their bulge/disk parameters indicate that these are disk-like sys- 
tems, with relatively weak bulges, making them unlikely candidates 
of major merger induced formation. Nevertheless, they do follow 
the same correlations, involving the Sersic index, as the other bright 
galaxies. This discrepancy needs to be explored with a larger sam- 
ple of bright, but disky lenticulars. 

4.1.4 Disk correlations 

There exist clear correlations between the two main disk parame- 
ters -central surface brightness of the disks (^d(O)) and disk scale- 
length (ra)- for both bright and faint lenticulars (see Figure 8) 
although they occupy different regions of the plot. Lenticulars with 
large disks have a lower central surface brightness, on average. A 
clear anti-correlation is seen between these two disk parameters for 
bright lenticulars but the scatter is large. The linear correlation co- 
efficient between /Xd(0) and log r<j is 0.43 with a significance of 
99.37 percent. Faint lenticulars show less scatter with the correla- 
tion coefficient of 0.70 at a significance level better than 99.99 %. 
Similar statistically significant correlation was reported by Khos- 
roshahi et al. (2000b) for early-type spirals and Mollenhoff & Heidt 
(2001) for late-type spirals for these disk parameters. In Figure 8, 
we also plot as dots the central surface brightness of the disks 
(Md(0)) against the disk scale length (r d ) of galaxies from the K 
band observations of Khosroshahi et al. (2000b) of early-type spi- 
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Figure 6. Sersic Index n as a function of unconvolved bulge central surface 
brightness. Symbols are as in Figure 2. 



Figure 8. Unconvolved disk central surface brightness fi^ (0) as a function 
of disk scale length. Early-type spirals (as dots) are overplotted for compar- 
ison. Other symbols are as in Figure 2. 
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Figure 7. Sersic Index n vs. bulge-to-total luminosity ratio. Symbols are as 
in Figure 2. 

rals. From the plot it is apparent that the disk scale lengths of our 
sample and the Khosroshahi et al. (2000b) sample span the same 
range, but pid(0) are brighter, on the average, for early-type spirals 
than faint lenticulars. 



4.2 Environmental dependence 

The dependence of cluster environment on formation and evolu- 
tion of galaxies is well known. The environment influence is best 
visible through the well-known morphology-density relationship, 
according to which early-type galaxies are preferentially found in 
rich environments (Dressier 1980), with spirals are found predom- 
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Figure 9. Dependence of the mean surface brightness within bulge effective 
radius ((^*b(< r c))) on bulge effective radius r e as a function of environ- 
ment. Faint field lenticulars are shown as open triangles while faint cluster 
lenticulars are plotted as filled squares. 



inately in low density environments. One of the stronger trends 
seen in the morphology-density relation at low redshift is the dra- 
matic increase in the number of lenticular galaxies in rich clus- 
ters (Dressier et al. 1997), suggesting that the dominant process 
in defining the morphology-density relation is the transformation 
of spiral galaxies into lenticular galaxies within rich clusters (Pog- 
gianti et al. 1999; Kodama & Smail 2001). These spiral galaxies 
are continuously supplied by accretion from the surrounding field 
during the course of the assembly of the cluster. Hubble Space Tele- 
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Figure 10. Unconvolved disk central surface brightness ^td(O) is plotted 
against disk scale length r<j as a function of environment. Symbols are as 
in Figure 9. 



scope (HST) imaging of 10 clusters at higher redshift (z ~ 0.3-0.5) 
by Dressier et al. 1997 shows similar behaviour. The mechanisms 
proposed for morphological transformation of spirals to lenticulars 
include dynamical interactions such as galaxy harassment (Moore 
et al. 1996, 1998, 1999) and gas dynamical processes e.g. ram pres- 
sure stripping (Abadi et al. 1999). Simulations have shown that 
galaxy harassment is a very efficient mechanism for transforming 
early-type disk galaxies to lenticulars (Moore 1999) with conse- 
quent reduction in disk sizes. Ram pressure stripping truncates the 
star formation by removing the cold neutral gas reservoir causing 
the disk component of a spiral galaxy to fade within the cluster 
environment. 

To address these possibilities observationally, it is necessary to 
investigate how scaling relations defined by various galaxy proper- 
ties vary between field and cluster environments. In this section, we 
examine a few correlations as a function of environment. Unfortu- 
nately, we are unable to include bright lenticulars in this compari- 
son because only 3 cluster galaxies out of the 38 in our sample, are 
bright. For faint galaxies, the number of galaxies in the field and in 
clusters are almost equal, making a comparison possible. 

The Kormendy relation for faint lenticulars plotted in Figure 2 
shows two distinct groups, unlike bright lenticulars. To investigate 
the origin of this dichotomy in the faint lenticulars, we plot the Ko- 
rmendy relation for faint lenticulars, differentiating them by envi- 
ronment, in Figure 9. We indicate field lenticulars with open trian- 
gles and cluster lenticulars with filled squares. The field lenticulars 
seem to follow the Kormendy relation reasonably well while the 
cluster lenticulars show a clear downward scatter with respect to the 
relation. On average, this requires a 1.5 mag arcsec~ fading of the 
mean surface brightness within r c of faint cluster lenticulars, with 
respect to the Kormendy relation. Boselli & Gavazzi (2006) have 
suggested that low-luminosity lenticulars in clusters might be the 
result of ram pressure stripping of late-type galaxies, which causes 
fading and thus lowers surface brightness. 

It is interesting to see if this fading is also seen for the disk 
component. In Figure 10, we plot the disk central surface bright- 
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Figure 11. Dependence of the r c — relation on the environment in K 
band for the faint lenticulars. Late-type spirals (as asterisks) are overplotted 
for comparison. Other symbols are as in Figure 9. The dotted line is the 
best fit to the faint lenticulars (field as well as cluster) excluding the one 
obvious outlier while the dashed line is the best fit to late-type spirals. Note 
that these two lines have similar slopes but a different intercept. 

ness against the disk scale length. As above, we indicate faint field 
and cluster lenticulars with different symbols. The faint field lentic- 
ulars exhibit a clear anti-correlation, whereas the cluster lenticulars 
occupy a limited region of the plot and show a downward scatter 
(indicating fading of the disk), very similar to the scatter seen in 
Figure 9. Our results obtained using only photometric data are con- 
sistent with the spectroscopic results of Barr et al. (2007), which 
support the theory that lenticular galaxies are formed when gas 
in normal spirals is removed (from both bulge and disk), possibly 
when well formed spirals fall into the cluster. It must be noted that 
the sample of Barr et al. (2007) is a subset of the BAM06 sample, 
which is included in our analysis. 

If lenticulars in clusters are indeed transformed spirals, it is 
likely that they preserve other signatures of their earlier existence. 
For instance, if they contained pseudobulges that formed through 
secular evolution, their bulge and disk sizes should be correlated 
(Courteau et al. 1996). In Figure 11, we plot the bulge effective ra- 
dius against the disk scale length for faint lenticulars (a similar plot 
for bright lenticulars may be found in Barway et al. (2007)). Faint 
field and cluster lenticulars seems to both follow the same correla- 
tion, but inhabit different regions of the plot. We also overplot data 
for a sample of late-type disk galaxies from Mollenhoff & Heidt 
(2001). The correlation for faint lenticulars and late type spirals 
has similar slopes (~ 1.40 ± 0.14 ) but different intercept which is 
expected as late-type spirals have larger disk scale length compared 
to faint lenticulars. This is a natural expectation if late-type spirals 
are transforming into faint lenticulars due to interaction with the 
cluster medium as well as with other galaxies in the cluster. 



5 SUMMARY 

The main results of this paper are: Several correlations such as the 
Kormendy relation, photometric plane etc. support the hypothesis 
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that bright and faint lenticulars (in cluster or field environments) are 
fundamentally different, with different formation histories. Bright 
lenticulars resemble ellipticals and bulges of early-type spirals sug- 
gesting that they may have formed like them - at early epochs 
via major mergers or rapid collapse. Faint lenticulars, on the other 
hand, have properties consistent with them having formed via in- 
ternal secular evolution processes (in the field) or via environment 
influenced secular evolution processes such as minor mergers, ram 
pressure stripping and galaxy harassment (in clusters). Although 
the dominant differentiating parameter between the two lenticular 
classes is luminosity, the environment also seems to play a role in 
determining the details of lenticular formation. In particular, the 
cluster environment seems to induce a fading of the bulge and 
the disk and possible transformation in morphology from spiral to 
lenticular, at least for the faint lenticular population. 

If the formation scenario of bright and faint lenticulars is in- 
deed completely different, it should also manifest as differences 
in star formation history. These differences can be probed using a 
combination of population synthesis models and multiband bulge 
disk decompositions. We intend to carry out such an analysis in a 
future project. 
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